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Abstract 

We associate some (old) convergent series related to definite inte¬ 
grals with the cyclotomic equation x m — 1 = 0, for several natural 
numbers to; for example, for m = 3, x 3 — 1 = (x — 1)(1 + x + x 2 ), leads 

t0 Jo dx {i+x+x*) = (T7J) = ( X - 3) + (! “ s) + (! “ I) + ■ ■ ■ • 111 
some cases, we express the results in terms of the Dirichlet characters. 
Generalizations for arbitrary to are well defined, but do imply integrals 
and/or series summations rather involved. 


1 Introduction 

Nicola Oresmes proved ca. 1350, the divergence of the harmonic series 
(e.g., see [1], p. 183). Indeed, today we know precisely how it does diverge 
([2], p. 14; see also [3]): 


(1.1) Harm(N) := — 

n=1 

For N » 1 

(1.2) HaxmfN) ^ + 7 + Y _ __1_ _ ... 
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where 7 is the Euler-Mascheroni constant (written and named C at times; 
for a recent reference, see [4]) 


(1.3) 7 = lim ( V - - log(N) ) « .577215 ... 

N ^°° V^I n / 

We do not know, even today (winter 2015), whether 7 is rational or not. 

Later, in 1668 Mercator ( = Kremer) proved ([1], p. 185) the convergence 
of the alternative series (of even/odd numbers) and summed it (i.e., taking 
the N —> 00 limit): 


(1.4) 


111 1 1 r 1 dx 

2 + 3 ~~ 4 + " ' + 21 V- 1 ~~ 2N + " ' y J 0 1 + x 


log{ 2 ) 


Apparently, some work form India preceded the (even later) so-called 
Gregory-Leibniz formula ([1] again, p. 184) for another alternative series (of 
inverse of odd numbers) 


(1 

In this paper we interpret (1.4) and (1.5) as arising from the cyclotomic 
equation of roots of the unity (Gauss; see, e.g. [5]): 


111 1 1 

1 ~3 + 5 _ 7 + "' + 41V-3~41V-l 


dx 


1 + x 2 


7 r 
4 


x m — 1 = 0 m € N; m = 2,3,4,...; 

(1.6) x m - 1 = (x - 1)(1 + x + x 2 + ... + x m-1 ) 

x 2m - 1 = ( x _ 1) (l + x 2 + x 4 + ... + x 2(m-l)^ 

In our cases, m = 2,4 respectively for (1.4) and (1.5), and then perhaps 
we can write anew the results in terms of some natural arithmetic functions. 
In (1.6), the roots / ±1 are complex conjugate pairs, as (1.6) describes a 
real equation, so the roots have modulus 1 , and all of them lie in the complex 
unit circle ~ S 1 . 

The main purpose of this paper is to generalize these results for other 
(generic) natural numbers m £ N. We shall try to relate the series to some 
definite integrals. 

We anticipate already here part of the workings for m = 2 and 4. 

For m = 2, x 2 — 1 = (x — l)(x + 1); we always leave out the x = 1 root. 
The inverse of (x + 1) enters in (1.4) and we proceed to the four operations 
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a) inversion , (x + 1 ) = (1 + x) —> (1 + x) 1 . 

b) expansion , (1 + x )^ 1 = 1 — x + x 2 — x 3 + x 4 — ... (|x| < 1 ). 

2 3 4 

c) term-by-term integration , —> x — ^\ + • • • 

d) taking the limits x = 1 minus x = 0 , —> (1 — |) + — |) + (| — ... 

Then we obtain the above result for m = 2: 

(1.7) 

(x + l) = (l + x)=> (l-x + x*-...)^ n (*- 4 ) + (t-X ) + - 

(x = 1 minus x = 0 ) = ^ + = log( 2 ) 

in terms of the arithmetic modulated “sign” function (=even/odd) /(n) := 

-—-—. Note that after the integration, the limit x = 1 can be taken. (1.7) 
involves the series 


(L8) E (^T ~h)= l dx TTx = log(2) 

expressed as a genuine definite integral. 

Similarly, for m = 4, it is x 4 — 1 = (x 2 — l)(x 2 + 1), separating the 
roots ±1 and ±i; now (1 + x 2 ) -1 enters in (1.5), which becomes, again after 
inversion, expansion, integration and taking limits: 


(1.9) 


(1 + x 2 ) 


2 \ —1 


/-I 1 \ /I 1 \ /I lx / 

( 3 ) + ( 5~7 ) + ( 9~IT ) + "' (= 4> 


which can be written as the (alternative) difference between the two series, 
and of course also as one integral 


(L10) E (4^3 “ sru) = l dx TTx* = I 

n=l x x 

And can be given now in terms of the so-called Dirichlet characters (see 
below). 

In this paper, as said, we shall generalize the constructions above (for 
rn = 2,4) for any integer rri £ N (in principle), and include (when possible) 
the appropriate Dirichlet character. 


3 



But first we want to fix our notation. For each m G N, m > 1, we shall 
consider first the finite sums, up to mN terms; so we define primarily 


(1.11) E m (JV) = l + - + . 


■ + ■ 


mN — (m — 1) mN — [m — 2) 


+ ...+ 


mN 


For example, £ 3 (N) = l + | + i + ... + -^ + ... + 3 ^_ 2 + 3^-1 + 3 lv> 
so these sums diverge , for N arbitrarily large, and indeed as (1.2) indicates: 


(1.12) Ti m {N) —» log(mN) +7 + -) = log(N) + 7 + log(m) 

m.JM 

(from now on we shall keep only the dominant and the constant terms in 
the divergent expansion ( 1 . 2 )). 

Also, we define partial bounded sums = Y? m (N) for % = 1,2,..., m) 
as displaced sums: 


(1.13) 


and 




1 


+ 


1 


m 1 ' m + 1 ' 2 m + l 
+ ...+ 


+ ...+ 


1 


(mN — (■m — 1)) 


y2 _ 1 | 1 

m 2 + m + 2 


(mlV — (m — 2)) 


(1.14) 


im-1 

-‘m 


1 


+ 


1 


m — 1 2 m. — 1 


+ ...+ 


1 


mlV — 1 


And finally, E”(-S + ^L. + ... + -J N ._ i-(l + \ + ... + jj) 


(1.15) (for IV »1) S-(AT)^^l + ^+o(l 

m m V A 


so 


Therefore, Y, m (N) (in 1.11) is the sum of all S^(A1) with i : {1, 2,... , m} : 


(1.16) S m (JV) = SumZA(N) = + ... + E£ 

Note these sums make sense, spite diverging, because all have only pos¬ 
itive terms and they are finite sums... 

The convergence or divergence of the above series is usually self-explained: 
convergence occurs always for the decreasing alternative series (Leibniz), but 
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as the convergence is not absolute, but conditional, the ordering in the se¬ 
ries should be maintained. The higher divergence behaves, if at all, with a 
factor oc log(N): divergences are no higher than logarithmic. We shall try 
to be careful in subtracting two divergent series... We use a philosophy close 
to physics (in particular, to quantum electrodynamics or q. e. d.): we first 
truncate the series, taking a fixed upper bound N 3> 1 (called the cut-off ; in 
physics this process is called regularization). Then we subtract other series, 
also divergent but with a similar type of divergence (this is called renor¬ 
malization): the result should be convergent (radiative corrections); see e.g. 
the book by J. Schwinger [ 6 ]. 

For a modern treatment of the Dirichlet series consult [7]. 


2 The case for m = 2 : Four methods 


Here we repeat the log( 2) result for the case m = 2 (1.7) by four different 
methods, because eventually the four might be useful. 

1) Redundance (direct) proof: ^) 2 (iV) and ^^(IV0> defined above, di¬ 
verge in a known way, for fixed (and large) N and, as said (1.11, 1.14) 

= S 2 - Si —► log(N) + log( 2) + 7 - hog(N) - 1 = 

( 2 . 1 ) 2 2 
= ^og(N) + ^ + log( 2) 

So 

( 2 . 2 ) —► (N —► oo) = log( 2 ) 

directly, without any integration or (infinite) series summation. 

2) Integration proof: Trivial here, as it is already (1.4): 

< M > fiTi = M2) 

3) Summation proof: we have the convergent series (so now N = oo) 
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which can be explicitly done using the digamma function 'I'(z) [see e.g. [8], 
pag. 258]: 


(2-5) 


* W :=«i^W)=m (^Z-:=Z\N) 


dz 


T(z) 


(for r(z) is the ordinary, Euler’s gamma function), with help of the expres¬ 
sion [8, p. 259, formula 6.3.16] 


OO 1 

— r~ -r = 4>(z) + 7 +- 

^ n(n + z ) z 

n =1 x 7 


for z = — \ . In total we get, as expected 


( 2 . 6 ) 


E 

n— 1 



4 log (2) 


4) Use of Hansen formula: This formula, again depending on the Diganuna 
function T(z), reads (Cfr. [9]; notice the sum starts in zero; x, y, z are some 
parameters) 


(2.7) 


E 

n.=0 


(■nx + y ) 2 — z 2 


1 

2 xz 




y + z 


- T 


y- z 


We put our n' above in Y ], ,, ) 1n as in' = n + 1), because n runs now 

from zero; hence this is equivalent to (2.7) with the parameter values x = 1, 
y = | and z = So our sum is (with 'h(l) = —7 and T(i) = — 7 —21og(2)) 


( 2 . 8 ) 


S = 2 


*(!)-*( 2 


41og(2) 


as it should be. 


Of the four methods, the most common and “easy” is the second (in¬ 
tegration); it can, in principle (that is, if the integrand is known and the 
integration feasible), always be used. The Mathematica Computer Programs 
give many integrals and double sums also directly. 

Please note in this (initial) rn = 2 case that what appears as result is 
the sign functionf(n ) := ^ 7 —; later, for some m > 2, (including 

m = 4) this will become a Dirichlet character: Here we have, as resume (to 
repeat) 
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(2.9) 


£* - S|(= lo 9 (2 )) = 2 


(- 1 ) 


n+1 


_ ^jodd _ j^even 


n =1 


n 


so NO Dirichlet character this time. 


3 The case for m = 3, 4, 6 


We discuss m = 3 first. Now (x 3 — 1) = (x — 1)(1 + x + x 2 ) as the 
roots are +1, oj = exp and conjugate Co; here u> is equivalent to a plane 
rotation by 120°. By direct integration, we obtain at once 

(3.1) [' to , = 

Jo 1 + X + x 2 3\/3 

which is a particular case of 


J — dx for A = (x — (jj){x — u>) = x 2 — (Troj)x + 1, Tru> = uj + oj 

For a general positive integer m we obtain, after an easy calculation (w 
is root of x m — 1 = 0) 


(3.2) 



dx 

1 — T rux + x 2 


it m — 2 
2 sin(^) m 


This formula (3.2) can be used at once for m = 3,4, and 6, yielding (3.1) 

and 


(3.3) 



dx 

1 — Tru ■ x + x 2 


— (m = 4) and =—— (m = 6). 
4 v 3>/3 


Later, we shall need also I = dx ,.j for cu = exp^^, for r 

natural between 1 and m — 1. The integral is now 


(3.4) 


7 r (m — 2r) 

( 2s in(^)) ' m 


(1 < r < (m 


1 )) 


Coming back to m = 3, the specific series is 1 — 5 + S' 2 — S 3 + S 4 — ..., 
with S = x + x 2 ; developing, we get (after term-by-term integration, and 
taking limits 1, 0) 
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1 


1 


(3.5) 


1 1 1 

^ + M 5^ + "' + 3iV-2 3JV-1 


E 

n =1 


i 


i 


3n — 2 3n — 1 


= v 1 _ v 2 — 

_ 2j 3 Zj 3 — 


7r 

3\/3 


Hence J [) 3 — ^ 3 finishes the calculation for m = 3, (as the third sum, 
X ) 3 trivially computable). 

( o\ 

In terms of the [10] Dirichlet character mod 3, namely x 2 j we have 


Ej - E? = -4= = X 

n=l 


00 (3)/ \ 

X 2 K n ) 


3y/3 


n 


where, of course XtEh 2,3,4, 5, 6 ) = (1, —1,0,1, —1,0) (so 3-periodic). We 
can even use the Dirichlet’s L-functions, which will include the j factors, 
but we refrain of doing that, as it does not illuminate the matter any further. 


As recapitulation, the series are obtained from the polynomial of roots 




(3.6) (m = 3) P 2 (x) = x 2 + x + 1 

after INVersion, EXPansion, INTegration and TAKing x = 1. In principle, 
the result of the series summation can be also obtained from the Hansen 
formula (2.7). 

Finally for this m = 3 case, we quote the four diverging series for com¬ 
pleteness: 


(3.7) 




1 


1 


72—1 


3n — 2 3n — 1' 


7r 

3\/3 


X 

? 7=1 


(3) 

X 2 


(n) 


n 


as the complete solution for the m = 3 case. We quote the following four 
divergent series for later use, still in this m = 3 case (the limit —> meaning 
just N » 1): 


(3.8) S 3 — 1 + - + ... + —■» log(N) + 7 + log( 3) 

oil 11 7 

(3.9) E? = - + - + +- > -log(N) + 1 

\ ) 3 3 6 3 N 3 J 3 
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1 


7r 


(3.10) SJ-1 + 1 + 1 + ... + —► jlog(N) + 1 + ^05(3) + 


(3.11) E? = - + - + ... +--- 

v ; 3 2 5 31V - 1 


l!09(JV) + 2 + li»9(3) - 


Now, for m = 4, we have first the natural cyclotomic expression 


(3.12) x 4 - 1 = (x 2 - l)(x 2 + 1) / 

(we write I because this is not the only possible factorization to be used). 
Repeating the steps as before for m = 3, our first final result is here 


(3.13) 


s1-e3 = E 


n= 1 


i 


i 


4n — 3 4n — 1 


dx 


x 2 + 1 


7r 


= arctan( 1 ) = — 


with \' 2 ™ (1, 2,3,4,5, 6 , 7, 8 ) = (1, 0, —1,0,1,0, —1,0,) is periodic mod 

4 (and restricted multiplicative). Note also why do we get i = 1 and 3 in 
(not 2!) mod 4: the expansion is for ~ 1 — x 2 + x 4 — x 6 etc., so it 
is with even powers only, so with only odd powers after integration! 


We are done with this, as £ 4 , £4 and £4 are automatically obtainable. 

The second factorization of x 4 — 1 is obtained by separating only the 
x = 1 root: 


(3.14) x 4 — 1 = (x — 1)(1 + x + x 2 + x 3 ) II 

As B := (1 + x + x 2 + .T 3 ) contains the x = — 1 root, one writes B = 
{x + 1)(1 + x 2 ), where the integral can be computed at once (indeed, it is 
indicated above). The final result for factorization II is 

\ l ° 9 ^ 2 ) + | 

- £ 2 computable from the above 
4 case by writing the analogous 


(3.15) 


dx 


1 + X + X 2 + x 3 


which is a redundant result, as defines £4 - 
calculation I. In fact, we end up this m = 
to (3.8) to (3.11): 
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(3.16) 


^ 4 ~\ + l + --- + W 


\logW + \ 


(3.17) vl = l + l + ...+ 


1 


41V-2 


9^2 — ~^1°9{N) + ^ + 9 ^ 0 5 ( 2 ) 


(3.18) 


X\ = l + l + ...+ 


1 


31V- 1 

-log(N) + 7 + -rlog(2) + — 


(3.19) S| —► \log{N) + 2 + llog( 2) - | 

The non-prime structure of m implies only a difference calculation, as it 
was also the case for m = 3, prime. We have c/>(3) = 2, </>(4) = 2. 

So the full solution for the m = 4 case has one redundancy (two fac¬ 
torizations), and it is done (solved) once a single calculation is made, e.g. 

fo dx TT^ = 7T / i - 

Now it is the turn of m = 6 . But here there are also several factorizations, 
as for the m = 4 case above. The simplest is (perhaps) 


(3.20) Firstly is x 6 - 1 = (z 3 - l)(.x 3 + 1) I 

As (ITx 3 ) = (1 + x)(l — x + x 2 ) : the integral is easy, with this factoring: 


(3.21) 


/•' dx 
o 1 T x 3 


-l°g( 2 ) + 


7r 

3\/3 


The operations INV, EXP, INT and TAK x = 1 applied to the polyno¬ 
mial P 3 (x) = 1 T x 3 yield the infinite but convergent sum 


(3.22) 


sj - Ej = £ 


n= 1 


1 


1 


6 n — 5 6 n — 2 


^Qff( 2 ) + vr 


3\/3 


Notice again the jump , now by three: is due to the cubic x 3n terms i 


m 


1 

l-\-x 3 * 
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All other factorizations are therefore redundant ; we just write them: 


(3.23) (x 6 — 1) = (x — 1)(1 + x + x 2 + x 3 + x 4 + x 5 ) =>- Eg — Eg 


II 


(3.24) 


(x b - 1) = (x - 1)(1 + x + x ) =7> Eg - Eg 


III 


Redundancy arises because from (3.22) one obtains Eg, as Eg is directly 
computable (our Sect. 2); so, in II, we know already the result, as Eg is 
again trivial; same in III, as Eg is again directly computable. 

It is remarkable here that in the original expansion (x 6 — 1) = 

(x 3 — l)(x 3 + 1) the simplest factorization (x 3 + 1) = (x + l)(l — x + x 2 ), 
yields for P 2 (x) = 1 — x + x 2 , the double series expansion 

(3.25) 

(1 + -) — (- + -) + (- + -) — (— + —) ± ... = Eg + Eg — Eg — Eg = 
v ^ 2> v 4 5 y v 7 8 ; v 10 11 6 6 6 6 


dx 


1 — x + x 2 


2?r 
3\/3 


which again it is a redundant calculation, as included in the (summed) 
Jo dx TTP- 

To finish these simple cases, we just write down the nontrivial summa¬ 
tions: 

1, . . . log( 2) login) logi 2) 7r 

S l(N) = -log(N) + 7/6 + + — 


(3.26) 


6 

4 _ log(N) 7 , M 3 ) 
6 + 6 + 4 


4\/3 


E 4 — 
^6 — 


7r 


12\/3 


l ^5 _ M3) 


y 1 _ yo _ 

Z^g Zjg — 


And, in this m = 6 case, we get also 
S 3 -S| = ilog(3) + ^. 
computing Eg and E| as in Sect. 1. 

In terms of some Dirichlet functions we get 


(3.27) 


S-S =E 


( 6 ) 

X2 

n 


, as (1,2,3,4,5,6) = (1,0,0,0,-1,0) = 


log(3) 
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For a modern treatment of sums involving harmonic numbers, see [4]. 


4 General numbers 


Here is the cyclotomic equation is, for m = 5, 

x 5 — 1 = (x — 1)(1 + x + x 2 + x 3 + x 4 ) = 

(4.1) _ 

= (x — l)(x — u)(x — u>)(x — UJ Z )(x — UJ 2 ) 

where u = exp = Rotation by 72° and co 2 = exp ^p = Rotation by 144°. 


We obtain easily 


(4.2) 


cos (72°) 


y/E- 1 

4 


cos (144°) 


-1-y/E 

4 


The integral T = f ( ] i + t . + j a (j3 +r 4 can be done as the real denominator 
splits in two real ones and quadratic, but we abstain to write explicitly 


(4.3) (1 + x 2 + x 3 + x 4 ) = (1 — xTru + x 2 )(l — xTrJ 2 + x 2 ) 

The expression for the integral is too long to write it. In summation 

terms > ^ is i+x+J+^i =■■■ = ( 1 -|)+(g-7) + - • • = £n=i (sTTi - 5iWs) 
still computable (e.g. with Mathematica), but still too long. The final result 
will be 


(4.4) 


dx 


-^- 5 - 7 = El - si 

1 + X + X 2 + X 3 + X 4 


This does NOT correspond to none mod 5 Dirichlet characters, but one 
can always put a (Dirichlet) function 


f~ 2 ^ (1, 2, 3,4, 5,6, 7...) = (1, —1,0, 0, 0,1, —1) periodic mod 5 
The simplest summation (according to Mathematica ) is 


(4.5) 


Vl v4 _ (! + \/5)tt 

5J2(5- >/5) 
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Notice the factorization in (4.3) does not allow us to write it as difference 
between two series, because the coefficients are not integer numbers. 

For the next prime, namely m = 7, we have three couples of complex 
roots, plus the x = 1 value: the sextic integral has not been attempted, but 
the summation can be again being done; we refrain to elaborate. 

This is the general trend for prime numbers p\ there are pairs 

of complex conjugate roots; “a priori”, the only integral vs. series is the 
simplest case, generalizing the above result 


(4.6) 



dx 


1 

1 + X + X 2 + . . . + xP~ L 



5 General case 


We deal now with some composite numbers. Compound numbers are 
easier; we just add a calculation for m = 8, with <f>(8) = 8 — 4 = 4: 


(5T) 

x 8 —1 = (x 4 —l)(x 4 +l), and (x 4 +l) = (x 2 —v / 2x+l)(x 2 +v / 2x+l), and so 


dx 


IT + log (3 + 2\/2) 


(5.2) 


x 4 + 1 


E 

n= 1 


4\/2 

1 


= . . . = So — So = 


1 


8n — 7 8n — 3 


7T log (3 + 2y/2) 


4\/2 


4y/2 


Another factorization is (x 8 — 1) = (x 2 — 1)(1 + x 2 + x 4 + x 6 ). It is 
equivalent to 


(5 ' 3) f 0 dX (l + x 2 + x 4 + x 6 ) ^(8n-7 8n-5) 

/ n=l x 7 

Still, a third factorization is (x 8 —1) = (x —l)(l + x+x 2 + .. .+x 7 ), which 
yields Eg — Eg, so the remaining Eg is obtained by difference (with Eg’ 4,6 ’ 8 
inmediate). The case m = 8 is potentially resolved. Again, we refrain to 
elaborate. 


13 



For m = 9, Eg and Eg can be computed directly, whereas the odd cases 
Eg’ 3,5,7 require further work, but it is again feasible. Similarly for m = 10. 

As general conclusion, we have showed a remarkable relation between 
the cyclotomic equation x m — 1 = 0 and some series and definite inte¬ 
grals; they go from the simplest integrals (and series) in the literature (like 
f dxj yy = log(2)) to very complicated cases, still feasible: the integrals 
have denominators factoring in quadratic ones, and the series are of the 
type —-j-L—_, computable, in principle, by means of the Hansen’s for¬ 
mula. 

There are, however, some questions left over in our work: for example, 
the series for 1 _. ) . + . i .j + 1 l ., {+ +rQ we identify with the series )T) 7 +1 ; this is cor¬ 
rect, but we have checked it “case by case”, offering no general proof; etc. 
Also we feel some new series might perhaps appear, whenever the quadratic 
components offer an integer expansion: those are two questions for the fu¬ 
ture... 
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